Mesoscopic thermoelectric transport near zero transmission energies 
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We study the thermoelectric transport coefficients of a one-dimensional (ID) electron waveguide 
connected to one and then two off-channel cavities, in presence of dephasing phonons. The model 
system is that of a linear chain described as a ID lattice. For simplicity we consider single-mode 
cavities, which may be tuned with external gates. While the presence of only one off-channel 
cavity yields a nearly symmetric transmission profile, which is cancelled around the cavity mode, 
an additional cavity modifies this profile strongly and results in an asymetric shape characterized 
by oscillations. In both cases, we consider electron-phonon interactions in our calculations and 
analyze their effects on the transmission function around the Fermi energy. Knowledge of the 
energy-dependent transmission function allows the numerical computation of thermoelectric trans- 
port coefficients, including the thermopower. In the presence of a second off-channel cavity, the 
sign of the thermopower depends on the relative position of this cavity energy level with respect 
to the Fermi energy: the thermopower is positive when low-energy electrons in the vicinity of the 
Fermi level are not transmitted, and becomes negative when the higher-energy electrons are not 
transmitted. 

PACS numbers: 72.10.-d, 72.15.Jf, 73.23.-b 



I. INTRODUCTION 

Thermoelectric systems are heat engines which, by us- 
ing the electron gas formed by their conduction electrons 
as a working fluid, directly convert a heat flux into elec- 
trical power and vice-versa, depending on the desired 
operation mode. In the linear regime, the coupling be- 
tween heat and electricity is embodied in a parameter 
called thermopower or Sccbcck coefficient, S. The prin- 
cipal virtues of thermoelectric heat engines are that their 
conversion efficiency is not size-dependent and that their 
operation does not rely on moving parts; the main draw- 
back is their modest energy conversion efficiency, which 
is at best of the order of 10 % of the efficiency of the ideal 
Carnot thermodynamic cycle. 

Since the first works reported by Seebeck [l[, Oersted 
0], and Peltier Q, general interest in thermoelectricity 
has much varied with time. Some significant milestones 
in this field are Thomson's thermodynamic analysis (if, 
Callcn's [U adaptation of Onsager's formalism Q, and 
Ioffe's works on device operation Q. In particular Ioffe 
related the thermoelectric conversion efficiency to a di- 
mensionless parameter denoted ZT, where T is the aver- 
age temperature across the device and Z is a parameter 
that depends on the material's properties, namely the 
thermopower S, and the electrical and thermal conduc- 
tivities, a and re respectively: ZT = aS 2 T / re (re entails 
both electron and lattice thermal conductivities, re e and 
Ki a t, but we neglect the latter throughout the present 
work since we are mainly interested in the electronic 
transport properties H). Though ZT has a profound 
meaning on the thermodynamic level, this parameter is 
much used as a means to measure how well in terms of 
energy conversion efficiency, a device operates. A good 
performance is associated to high values of ZT (typically 



around 1). 

Enhancement of the performance can be obtained fol- 
lowing three routes: i j optimization of the thermoelectric 
material's properties H; iij optimization of device work- 
ing conditions through appropriate design fs, [ljj and 
electrical and thermal impedance matching s 146 1; a nd Hi/ 
design of new nano- and mesostructures [12L Il3| . The 
first route, which essentially consists of finding means to 
decrease re and increase a has produced a number of very 
interesting results, but further significant progress seems 
out of reach 14|, one of the difficulties being that both 



a and re are tightly linked; this link is the Wiedcmann- 
Franz law for metals. The second line of work is based 
on the observation that the highest performances are usu- 
ally obtained in a very limited temperature range for a 
given material, so that devices must be designed in such 
a way that they may operate over a large temperature 
range, which is usually achieved by segmentation; in ad- 
dition, when seeking maximum efficiency at maximum 
output power for a thermoelectric generator, one should 
ensure that the conditions for both thermal and electrical 
impedance matchings are fulfilled. 

The third route, which we explore in this paper, was 
permitted thanks to the rapid and tremendous progress 
in the field of fabrication of nano- and mesoscopic ar- 
tificial structures. The interest of going down to such 
scales stems from a number of effects that truly enrich the 
physics of thermoelectricity and pave the way to possi- 
bly better device performance; we can mention quantum 
confinement and the related modification of the carriers' 
density of states, mesoscopic fluctuations, breakdown of 
the Fermi liquid picture, and limits of validity of standard 
thermodynamic approaches due to finite-size effects. Fur- 
ther, since for macroscopic systems an increase of ther- 
mopower impacts on the values of thermal and electrical 
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conductivities and may not necessarily result in a signif- 
icant increase of ZT, it is also interesting to see how at 
the mesoscopic level the quantum effects pertaining to 
electron confinement and quantization of the transport 
coefficients [HI, [l6| may modify the interplay between S, 
ft, and a. In this paper we see how to achieve high figures 
of merit with systems having large thermopowcr and low 
thermal conductivity at the same time. 

We focuse on a one-dimensional (ID) system in the 
presence of off-channel cavities, which perturbe the 
electrical conductivity a or, equivalently, the energy- 
dependent transmission function T. This system 
may be fabricated from structures which contain two- 
dimensional electron gases. The main property of the 
mesoscopic system we study is the cancellation of the 
transmission function at a given energy we may choose 
and vary using an external gate; other characteristics are 
described in Section 2. Now, the main idea underlying 
the present work comes from the Cutler-Mott formula 
[l7l | , which gives the Seebeck coefficient at low tempera- 
ture: 



S = - 



7T 2 k%T fd\n(a 



dE 



(1) 



E=E F 



where the derivative of the electrical conductivity is taken 
at the Fermi energy E-p. Equation ([TJ which is valid for 
metallic systems such as that under consideration in this 
paper, shows that when the transmission vanishes, the 
Seebeck coefficient tends to infinity. In general, when 
the transmission vanishes, the derivative da /dE vanishes 
too, but in virtue of l'Hopital's rule, we know that the 
limit does diverge. Therefore, if the conduction band 
contains a transmission-cancelling energy level, the cou- 
pling between heat and electrical fluxes is enhanced in 
the vicinity of this energy, which in turn may yield a 
significant increase of the figure of merit ZT. 

As will become obvious with the definitions of the 
transport coefficients given further below, a symmetric 
transmission function would tend to cancel the contribu- 
tions from higher and lower energy terms with respect 
to the Fermi level, and hinder transport, so we seek con- 
figurations that makes the transmission function profile 
asymmetric. Indeed, the transmission function should 
be asymmetric around the Fermi energy so that the hot- 
reservoir electrons may leave their high energy states and 
proceed against the applied source-drain bias through 
the waveguide to the cold reservoir's available states as 
noted in Ref. [l8| . In this work, we investigate the effects 
of electron-phonon coupling on the transmission profile, 
and see whether this coupling impacts positively on the 
transport coefficients and figure of merit. 

Our phenomenological model is based on a number of 
assumptions. First, our main framework is that of single- 
particle theory so that the ballistic electron transport in 
the waveguide may be described by a single-band effec- 
tive mass Hamiltonian. We consider the linear response 
only, which implies that the temperature gradient across 




FIG. 1. (color online) An off-channel cavity connected to ID 
electron waveguide. The shape of the cavity and its coupling 
to the waveguide are controlled by external gates. 



the system is small; incidently this assumption provides 
a small energy scale for the model. The temperature T 
is sufficiently small so that the Cutler-Mott formula (JTJ 
remains valid. Effects of the couplings of the waveguide 
to the off-channel cavities and the phonons are studied 
with the nonequilibrium Green's function approach fl9| . 
adapting to our purposes recent works on dephasing ef- 
fects on quantum transport (2114241 ] . We assume that the 
systems operates in a steady-state regime, and that the 
two reservoirs, to which the electron waveguide is con- 
nected, are perfect. The Fermi energy E-p in the whole 
system is given by the highest level occupied by an elec- 
tron in either reservoirs. The phonons are assumed to be 
in equilibrium. 

The article is organized as follows. In section II, we 
present in detail the basic system under consideration. 
The system Hamiltonian, the thermoelectric transport 
coefficients and the transmission function are defined and 
discussed considering the presence of only one off-channel 
cavity coupled to one site of the lattice. In section III, we 
introduce the effects of the coupling of parts of the ID 
lattice to a bath of dephasing phonons; for this, we rede- 
fine our basic chain model using a decimation procedure. 
In section IV, we study the response of the system's ther- 
moelectric transport coefficients to the presence of an ad- 
ditional off-channel cavity. The article ends on conclusive 
remarks, followed by the Appendix where some detail on 
the calculation of the self-energies using the decimation 
procedure is given. 



II. MODEL OF THE BASIC SYSTEM 

We consider a ID electron waveguide, which we de- 
scribe as a ID lattice; electron conduction is character- 
ized by hopping between the lattice sites. The waveguide 
is connected first to one off-channel cavity by coupling to 
one of the lattice sites. We also assume a small cav- 
ity size so that it contains only one energy level accessi- 
ble to the conduction electrons [25| and, in the absence 
of the off-channel cavity, the system is uniform and has 
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full transmission: T(E) = 1, \/E, at any Fermi energy. 
The cavity is controlled by external gates which handle 
its coupling to the electron waveguide and its shape as 
shown in Fig. [TJ 



A. Hamiltonian 

The Hamiltonian of the system can be written in a 
tight binding model as the sum of three contributions: 



H — Hf + Ha 



H c 



(2) 



where Hi is the electron waveguide Hamiltonian, which 
describes the transport of an electron from — oo to +oo 
in the absence of the off-channel cavity: 



H, = 



J2 -t(|x)<a! + l| + |a; + l)<x|)+co|0)(0| (3) 



where \x) is a Wannier state at location x along the ID 
lattice, and t is the hopping coefficient, which we take 
as: t = ti 2 /2m*a 2 in the effective mass approximation 
(m* is the electron effective mass, and a the lattice spac- 
ing). The energy eo is that of the site at the center of 
lattice (x = 0). The dispersion relation derived from the 
Hamiltonian Hg reads E = — 2icos(fc), leading to a con- 
duction band E E [—2t, +2t]. The term H d characterizes 
the off-channel cavity with a single level: 



H A = V d \d)(d\ 



(4) 



where the energy V d of the state \d) can be controlled 
by an external gate. The term H c characterizes the cou- 
pling between the cavity and the site at the center of the 
electron waveguide to which it is coupled: 



H c = t c \0){d\+t c \d)(0\ 



(5) 



where the coupling parameter is noted t c . 
This simple model may represent a multilevel quantum 
dot when the level spacing is large compared to the range 
of energies contributing to the transport due to tempera- 
ture smearing. In all that follows, the hopping coefficient 
in the leads is taken as t = 1 such that all the energies in 
our study are expressed in the units of t. 



B. Thermoelectric transport 

We consider three transport coefficients: the electrical 
conductivity <r, the thermopower S, and the electronic 
thermal conductivity Kq, which is related to the zero- 
electric-current thermal conductivity k c via: 



kq = k c + TaS 



(6) 



These three coefficients are given by the standard expres- 
sions of the solutions of the Boltzmann equation [261 ] : 



TaS = f{E-n) f-^jT(E)dE 



Tk 



3/ 
dE 



T(E)dE 



(7) 



(8) 



(9) 



where /i is the electrochemical potential of the electron 
system. The important quantity that determines all the 
necessary physical quantities for our work is the trans- 
mission function T. Knowledge of the three transport co- 
efficients thus permits an analysis of the figure of merit, 
which may be rewritten as: 



ZT = 



aTS 2 1 Ko 
1 - cjTS 2 /k 



(10) 



Equations (J7J, (JSJ, and ((9]) show that the electron 
waveguide, via its transmission function T which we de- 
fine below, acts as an energy-selective filter [l^]. As 
mentionned in the Introduction, for electron transport 
to take place, it is necessary that the transmission func- 
tion be asymmetric, and this is clear with Eq. ((5J) : since 
(E — fx) changes sign around the Fermi energy, a symmet- 
ric transmission function would yield cancellation. Also 
note that (E — fj) corresponds to the heat involved in the 
transport process, so a high energy conversion efficiency 
may be obtained on the conditions that the transmis- 
sion profile permits a high elctrical flux from hot to cold 
reservoirs with a minimal heat flux. 



C. Transmission function 

The energy- dependent transmission function may be 
obtained from the Fisher-Lee formula 123 : 



T(E) = r L Gr R Gt 



(ii) 



where G is the retarded single-particle Green function of 
the central site (Wannier state |0) with energy eo), taking 
into account the effect of the leads and the off-channel 
cavity by means of their self energies: 



G 



1 



E 



E-V d 



(12) 



where the self energies of the left (El) and the right(ER,) 
leads are taken equal: Er = El = E, and have the 
following energy-dependence [28[: 



E 
~2 



(13) 



By "central" region, we mean the site |0) connected to 
the cavity and the leads at the same time. In order to 
have full transparency (T = 1) at all the Fermi energies, 
the energy eo should be equal to that of the lead sites: 
eq = 0. Note that we kept eo explicit in Eq. (|12p despite 
its vanishing value, in order to better see the partition of 
the system and define the left and right self energies. The 
quantities IY and Tr in (jllj) characterize the coupling 
of the central region to the leads, and they are seen as 
linewidths. Their expressions are given by the imaginary 
part of the self energy of the leads: 




1.7 



1.6 



FIG. 2. (color online) Transmission profile in the presence of 
the off-channel cavity (orange curve) and without the cavity 
(dashed blue). The presence of the cavity changes drastically 
the transmission of the system and can even lead to cancel- 
lation when the Fermi energy equals the cavity level (here 
V A = -1.8). 



L,R 



-23E 



L,R 



(14) 



This partition into a central system and the surround- 
ing leads and cavity is optional and not unique (29j . Our 
choice is so to obtain a scalar Green's function, which 
greatly simplifies the expression of the transmission func- 
tion. Notice in the expression of the Green's function 
(fT2"j). that the off-channel cavity is taken into account 
through its effects on the central site, which can be un- 
derstood from the decimation procedures [30|, HH ; these 
effects may be contained into a self-energy: 



ti 



E-V d 



(15) 



Substitution of the above quantities into the Fisher- 
Lee formula (jll[) yields the following simple expression 
for the transmission function of the system: 



T(E) 



l + [tl/T(E-V d )Y 



(16) 



The profile of the transmission is given in Fig. [2] We 
notice that the presence of the cavity changes drastically 
the profile of the transmission and it even leads to its 
cancellation, which can be viewed as a quantum destruc- 
tive interference. The value at which this cancellation 
appears is directly obtained from the transmission func- 
tion, in Eq. (fTtjj) : as E ->■ Vd, T — > 0, and the width of 
the antiresonance is controlled by the coupling parame- 
ter t c . A first-order expansion near the energy E = V d 
shows that the transmission function is parabolic: 



r(2S)~{~) (E-V A f 



(17) 



and it becomes clear that the weaker the coupling to the 
cavity is (t\ <C T ) the sharper the antiresonance is. To 
recover the full transmission of the system (as without the 
cavity) one may increase the potential Vd to infinity, since 
this implies that the cavity region becomes forbidden to 
the electrons moving across the waveguide. Decreasing 
the coupling t c down to makes the antiresonance nar- 
rower but does not eliminate it; and setting t c = is a 
rather abrupt way to take the cavity off. 

We just saw the effects of the coupling parameter t c 
and the cavity level Vd on the transmission function, and 
clearly they have no influence on the symmetry of the 
profile. In order to induce some asymmetry we need to 
introduce an additional physical process, hence a new pa- 
rameter. In what follows, we allow the electrons that flow 
through the waveguide to interact with phonons within 
two small regions located on both sides of the site coupled 
to the cavity. 



III. EFFECTS OF DEPHASING PHONONS 

In this section, to study the effects of electron-phonon 
interaction on the transport coefficients, we consider the 
presence of external phonon baths symmetrically located 
on both sides of the cavity and connected to n — 2 x N 
sites. The system remains essentially the same as that of 
the previous section where we considered phase-coherent 
transport, but we now assume that in some parts of the 
left and right leads, the electrons are interacting with 
phonons as shown in Fig. [3J The description of the new 
situation is done using a decimation procedure, which is 
described further below. 
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a) 



Connected to phonons bath 



Vd 



Connected to phonons bath 



b) 



-o— 



-o- 



Renormalized sites after 
decimation 



FIG. 3. (color online) Sketch of the decimation procedure (a) — > (b). 



A. Electron-phonon interaction 

Following Ref. [l9j . we make approximations to per- 
form tractable calculations: we assume that the phonons, 
which form a bath of independent oscillators, are in equi- 
librium; we consider only one-phonon scattering pro- 
cesses (self-consistent Born approximation) ; the electron- 
phonon interaction is local and has no effect on the hop- 
ping coefficient 20]. Further, we simplify our analysis by 
considering single- mode phonons, with energy Huiq. The 
phonon-electron interaction Hamiltonian thus assumes a 
simple form and reads: 



ph = A^|x)(o;|(6 + 6 t ) 



where and b respectively are the second-quantized 
phonon creation and annihilation operators, and A is the 
electron-phonon interaction strength considered here as 
site- independent . 

The approach to access the transmission function 
through the device in presence of dephasing phonons 
starts with expressing the retarded (advanced) and lesser 
(greater) self energies as done in Refs. [2lL [23l l24l. |32|. [33j . 
These quantities are obtained in a closed form within 
the self-consistent Born approximation [2l], [34[, which 
ensures current conservation [2lL ■ 



B. Decimation procedure 

To simplify the description of the system connected to 
the phonon baths on both sides of the central site, which 
is connected to the off-channel cavity, we use a decima- 
tion procedure [3(J [H| as depicted in Figs.[3Ja) and[^6): 
the red sites are connected to the phonon baths, and the 
green site is at the energy level of the off-channel cavity; 
the system thus is formed of two parts: the interacting re- 
gion connected on its left and right sides to semi-infinite 
perfect ID leads. The decimation consists in taking out 



the central part (five grey sites and the green one), and 
rcnormalizing the sites (and their links) to which it was 
attached. The renormalized sites are shown in blue. This 
procedure is exact in the sense that both systems (a) and 
(b) are equivalent for transmission since they share the 
same Green's function submatrix related to the remain- 
ing sites. The sites are renormalized as follows: we note 
h the Hamiltonian of the part we take out and we define 
the Green's function g = E 1 _ h ; then the Hamiltonian of 
the sites directly coupled to the removed part is: 



h, = h[V + rt 



1 



E-h 



(19) 



(18) where h T is the new Hamiltonian of the renormalized sites 
(a sub-Hamiltonian of the total one), hr is the Hamil- 
tonian of these sites in the absence of coupling to the 
central part, and r is the coupling matrix between the 
central part and its surrounding sites (t' is its complex 
conjugate). Equation (|T9| clearly shows that the Hamil- 
tonian hi is energy-dependent, and we recognize the form 



of a self-energy: r 



t. 



E-h 



Actually the decimation proce- 



dure may be seen as the use of the concept of self-energy 
in the formalism of Landaucr-Buttikkcr to take into ac- 
count the effects of the reservoirs, or infinite leads (an 
application of such procedure may also be in Ref. |36(). 
Here, since the part of the system which is taken out 
is finite, the self-energy is real valued. The aim of the 
decimation procedure is to lower the dimension of the 
Green's matrix since it is numerically obtained using a 
recursive algorithm. Another interest of this procedure 
is that it simplifies the equations describing the phonon 
self-energy E ph , entering the definition of the retarded 
Green's function: 



G 



I 



(20) 



E - H - S loads - £P h 

where, after application of the decimation procedure, 
the Hamiltonian H takes the form of an n x n ma- 
trix. The leads are characterized by their self-energy 
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matrix £ loads , which has only two nonzero elements: 
pleads = pleads = s _ The phonon se lf-energy is more 

complicated and it is obtained self-consistcntly with the 
Green's function. Within the Born approximation and 
the assumption that the effect of the phonons is local 
(their effect is restricted on single sites and does not 
change the hopping coefficients) the phonon self-energy 
reads: 



-'ph 



7 2 PG 



(21) 



where T> is a superopcrator whose application to a matrix 
sets all its elements, except the diagonal ones, to zero, 
thus returning a diagonal matrix 21|. The parameter 
7 characterizes the coupling between the phonon baths 
and the system, hence the dephasing process. Note that 
the Green's matrix ([2U| and the phonon self-energy (f2"Tj) 
must be computed using a recursive process starting from 
an initial value of E ph , the computation being continued 
until convergence is reached to a good precision. 



C. Transmission function 

As explained in the works of Cresti and co-workers [2lL 
HH H3 j the transmission T at energy E may be written 
as the sum of two contributions, which are called the 
coherent and incoherent terms: 



T(E) = Tcoh(E) + T mc (E) 



(22) 



The coherent term originates in the leads' lesser self en- 
ergies [2l| and assumes the following form: 



electrons through the cavity. The numerical simulations 
based on an iterative procedure yield the results plotted 
in Fig. 2J We discuss first the case without connection 
to the off-channel cavity. Comparison of Figs. [2] and 0] 
shows that the presence of phonons drastically perturbes 
the transmission function: the amplitude is significantly 
reduced and oscillations appear. An increase of the size 
I of the interacting regions reduces the transmission am- 
plitude of the system, which reflects the dephasing effect; 
an increase of the distance L between the two interacting 
regions induces more oscillations in the transmission pro- 
file: the system mimicks the behavior of a Fabry-Perot 
interferometer. Indeed, the presence of phonons renor- 
malizes the potentials of the lattice sites with which they 
interact, via the real part of the self energy E ph , and 
since there are two distant electron-phonon interacting 
regions, the system shows oscillations of its transmission 
function. In presence of an off-channel cavity, the trans- 
mission profile remains essentially the same, except near 
the cavity level where the transmission is cancelled. Note 
that the transmission always cancels at the cavity level: 
the electron-phonon coupling, as described in our model, 
does not impact on the location of the antiresonancc. 
The very interesting result we observe is the asymmet- 
ric aspect that the transmission profile acquires around 
the Fermi energy in presence of dephasing phonons. The 
transmission of either high or low energy electrons can be 
reduced depending on whether the cancellation occurs at 
the local maximum or minimum of an oscillation: the 
system then acts as an energy-dependent carrier filter. 



D. Thermoelectric coefficients in presence of 
phonons 



(23) 



and the incoherent term, in the phonon's lesser self- 
energies f2lj . with a similar form: 



Tinc{E) = Fl-BihTr 



(24) 



In Eqs. (|2"3")) and (pM)) . the couplings Tl and Tr are those 
defined in Eq. (fl"4")l . and the matrix elements Q\ n and B\ n 
(in the notations of Cresti [2l|) are defined as follows. 
The matrix elements of Q are: 



Qij — | | 

and the matrix B reads: 



B 



Q' 2 



1 - 7 2 



(25) 



(26) 



We may now analyze the effects of the coupling of the 
electron waveguide to phonons on the transmission of 



Once the profile of the transmission in known, the 
transport coefficients may be directly obtained by nu- 
merical integration of Eqs. (JT]), (O, and (J9j> ; these coeffi- 
cients are shown as a functions of temperature in Fig. [5] 
We choose to analyze two different situations: one cor- 
responds to a system more transparent to high-energy 
electrons, the other is more transparent to low-energy 
electrons. In both cases we notice that the conductance 
a increases with temperature (note that in the linear 
regime, the temperature variation must remain small). 
This can be understood as follows: since the antireson- 
nace at the Fermi energy is extremely narrow, a temper- 
ature increase allows electrons with energy around the 
Fermi energy within a fceT range, to be transmitted. 
The behavior of the electronic thermal conductivity K e 
is similar for the same reasons. 

Both profiles of thermopower S as function of temper- 
ature T show an extremum. The difference of signs is 
due to the location of the antiresonance of the transmis- 
sion function, which corresponds either to a minimum or 
a maximum of the oscillations around the Fermi energy: 
in the upper row of Fig. [SJ the system is more trans- 
parent to electrons with energy E < E-p whereas, in the 
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FIG. 5. (color online) Temperature dependence of the transport coefficients for two different distances between the electron- 
phonon interacting regions. Top panels: L = 26a; bottom panels: L — 16a; and for both cases I = 6a. The related transmission 
profiles are shown in FigfJ] 



lower row, it is more transparent to electrons with energy 
E > E-p. This change of sign for the thermopower could 
also be anticipated from the formula © . We remind that 
this asymmetry is due to the modulation created by the 
two electron-phonon interacting regions. Note that while 
the values of the figure of merit ZT are modest in both 



cases, their shapes indicate that with appropriate adjust- 
ments of the location energy (which corresponds to the 
peak) and passband, very interesting developments are 
possible since the energy spectrum of the thermoelectric 
system is tunable. Indeed, this is a key point as regards 
the thermodynamic coupling of the system to its environ- 
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ment or other systems: minimization of the mismatch of 
the respecyive energy spectra of two systems permits a 
lowering of entropy production as these two systems get 
coupled (this is well known in the case of monochromatic 
photons interacting with a black body (37|). 



IV. SYSTEM CONNECTED TO TWO 
OFF-CHANNEL CAVITIES 

The analysis we did so far concerned various situa- 
tions for which the electron-phonon interaction-induced 
asymmetry is at most typically of the order of the ampli- 
tude of the oscillations. We now see whether a stronger 
asymmetry is possible by simply connecting an additional 
off-channel cavity to the electron waveguide, still in the 
presence of dephasing phonons as in the previous case. 

The two cavities arc supposed to be small and each one 
contains only one energy level that is tunable with exter- 
nal gates. The presence of a second off-channel cavity 
adds a new parameter to the problem. When connected 
to the ID system, the additional cavity also induces a 
cancellation of the transmission at an energy equal to 
that of its level. We now assume that the Fermi energy 
is equal to the energy level of one of the cavities and let 
the other level free to be changed. The Hamiltonian of 
the new system is quite similar to that of the previous 
case: H = H e + H c + H d + H c _ ph , with 



+00 



H, 



-t(\x)(x + l\ + \x + l)(x\) (27) 



X — — OO 



H d = V d \d){d\ + e d \d'){d'\ 
H c = t c \0)(d\+t c \l){d'\+h.c 



(28) 
(29) 



where Hi is the Hamiltonian of the electron waveguide, 
H d is the Hamiltonian of the two cavities with V d and 
e d their levels, H c is the Hamiltonian that characterizes 
the coupling between the waveguide and the two cavities. 
For simplicity, we choose the same coupling parameter t c 
for both cavities. Here I represents the index of the site 
to which the second level is connected (the first cavity is 
connected to the site with index 0). 

The transmission function is obtained from Eqs. (f2"2"j). 
(|2"5*)) , and (f2"B")) . By varying the level e d , we change the 
profile of the transmission T and obtain situations where 
the low or high energy electrons in the vicinity of the 
Fermi level cannot be transmitted due to the overlap of 
the two antircsonances as can be seen in Fig. [5] The 
simulation is performed considering system with a non- 
interacting central region of n — 67 sites at a Fermi en- 
ergy E-p = —1.8 (equal to the cavity level V d ). 

In order to reach significant values of ZT, we investi- 
gate all the possible situations corresponding to various 
differences between the two cancellation energies, e d — V d , 
and we plot the figure of merit ZT as a function of the 
temperature T . The numerical simulations are based on 
the same procedure as for the previous section. We only 



T(E) 




FIG. 6. (color online) Transmission profile in the presence of 
two cavities. The level of the first cavity is chosen to be Vd = 
Ef = —1.8. The energy level of the other cavity is ed such 
that: (e d --E F )/to = ^ 2 - 24 (blue curve); (e d -E F )/t 2 c = -0.64 
(green curve); and (ej — E-p)/t 2 = 0.32 (red curve); all with 
t c — 0.025. The overlap of the two cancellations yields a large 
zero-transmission region. The green curve corresponds to the 
maximum of ZT in Fig(7] The central noninteracting region 
contains n = 67 sites and 7 = 0.81. The distance between the 
two cavities is 50a. 



have to change the Hamiltonian of the central system 
by the new one including the additional cavity and iter- 
ate the algorithm in order to obtain in a self-consistent 
fashion the Green's matrix G and the phonon self-energy 
£P h . The results are shown in Fig. [7J The figure of 
merit ZT has a maximum at ZT = 5.38 obtained for 
[e d —E-p)/t^. = —0.64. This maximum and its location de- 
pend on the system's characteristics (central region and 
the phonon's regions). We notice the presence of another 
local (but lower) maximum in the ZT curve. 

We also computed the thermopower 5 as a function 
of the same parameters (ed — E-p and T). The Seebeck 
coefficient has a maximum which corresponds to the max- 
imum of ZT and a minimum which corresponds to the 
local maximum of ZT factor. The sign of the Seebeck co- 
efficient changes as e d varies with respect to the Fermi en- 
ergy: the thermopower is positive when low-energy elec- 
trons in the vicinity of the Fermi level are not transmit- 
ted (see Fig. [7]) whereas it becomes negative when high 
energy-electrons are not transmitted. Note that it is also 
possible that for a same configuration of the system, the 
Seebeck also changes sign because of temperature smear- 
ing which tends to give more weight either to lower- or 
higher-energy electrons as the temperature increases. 



V. SUMMARY AND CONCLUSIVE REMARKS 

Our study of the thermoelectric transport coefficients 
of a ID electron waveguide connected to one and then two 
off-channel cavities, in presence of dephasing phonons, 
is based on the Landauer-Biittikcr formalism and the 
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FIG. 7. (color online) Transport coefficients as functions of the energy difference — Ef, and temperature T. 



noncquilibrium Green's function techniques. We showed 
how to obtain a strongly asymmetric shape for the trans- 
mission function: first by coupling parts of the waveguide 
connected to an off-channel cavity to phonons, and then 
by coupling this system to an additional off-channel cav- 
ity. In our case, the transmission profile is character- 
ized by oscillations, which trigger changes of sign of the 
thcrmopower as the location of the antiresonance of the 
transmission function corresponds either to a minimum 
or a maximum of these oscillations around the Fermi en- 
ergy. This sign change shows that the system acts as an 
energy-selective filter: low-energy electrons in the vicin- 
ity of the Fermi level are not transmitted (case of pos- 
itive thermopower) , and the higher-energy electrons are 
not transmitted (case of negative thermopower) . 

As pointed out in the Introduction, from a thermo- 
dynamic viewpoint a thermoelectric system is a thermal 
engine to which a ZT-dcpendcnt energy conversion effi- 
ciency is associated. Recently some aspects of the ques- 
tion of efficiency at maximum power of low-dimensional 
thermoelectric systems has been discussed in terms of 
performance for quantum dots and ID ballistic conduc- 
tors (38M401 ]. These works pertain to the more general 
framework of finite-time thermodynamics (4ll-|44| , a field 
which contains a number of interesting open questions 
at the macroscopic level pBI - l48j . which must also be ad- 
dressed at the mesoscopic level; these questions are re- 
lated to the location of irreversibility sources in the sys- 
tem as a whole and the impact on the efficiency at max- 
imum power, two particular cases being those of endore- 
versiblc and exoreversible engines. The question of irre- 
versibility at the mesoscopic level is certainly not trivial 
since this requires a careful characterization of the cou- 
pling of a system to its environment. 
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GREEN'S MATRIX DERIVED FROM THE 
DECIMATION PROCEDURE 

As an example, we give the matrices defining the 
Green's function for the model presented in Fig. [31 First 
we define the Hamiltonian h of the part we want to take 
out: 



(30) 



where the diagonal elements ha represent the potentials 
on the sites of the central part including the cavity level, 
and the off-diagonal elements h%j show if the site i is 
connected to site j and with which strength. 

To perform the decimation procedure, we take out this 
part and renormalize the sites to which it was attached. 
We first express the Green's function g = -grzr- After 
decimation, we obtain the system shown in Fig. [3t> with 
the Hamiltonian of the region between the leads reads: 



-1 0' 

-10-10 

-1 .911 gin 

g n i g nn -1 

-1 -1 

-1 



H 



(31) 



Note that we used the fact that the coupling of the central 
system to its environment (lead with phonons) is t = 1; 
otherwise t 2 gij must be substituted to gij. 
The self energy matrix of the leads reads: 
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where £ is the self-energy of each semi-infinite perfect 
lead given in Eq. (| 13[) . The recursion process that drives 



the numerical computation starts with the initial value 
of E ph = S loads and continues with equations (f^Tj) and 

(20). 
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